Constraining CP violation in neutral meson mixing with theory input 
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There has been a lot of recent interest in the experimental hints of CP violation in SjJ j, mixing, 
which would be a clear signal of beyond the standard model physics (with higher significance). We 
derive a new relation for the mixing parameters, which allows clearer interpretation of the data in 
models in which new physics enters in M12 and/or ri2. Our results imply that the central value of 
the D0 measurement of the semileptonic CP asymmetry in BjJ^ decay is not only in conflict with 
the standard model, but in a stronger tension with data on ATs than previously appreciated. This 
result can be used to improve the constraint on AT or Asl, whichever is less precisely measured. 



I. INTRODUCTION 

Recently, CP violation in neutral meson mixing re- 
ceived renewed attention due to the D0 hint of CP vio- 
lation in B-B mixing, measured by the CP asymmetry 
in decays of a hh pair to two same-sign muons m , 

^SL = -[7-87 ± 1.72 (stat) ± 0.93 (syst)] x 10"^ . (1) 

At the Tevatron both B^ and i?j are produced, and hence 
AgL is a linear combination of the two asymmetries p| 






(0.594 ± 0.022) A^L + (0-406 ± 0.022) A^l • (2) 

The central value in Eq. ([Ij would be a clear sign of 
new physics (NP) [2|, S]- Measurements at the e+e" B 
factories [1] and at D0 [i] yield 



Ad 



-(0.5±5.6)xl0" 



A: 



SL 



(1.7±9.2)xl0"3. 

(3) 

In i?^ mixing, a nonzero lifetime difference, AF^ — 
Tl — Th , was established recently. 



AF, 

AF, 
AF, 



(0.116 ± 0.019) ps" 
(0.068 ± 0.027) ps' 
(0.163l°:g65)ps-i 



LHCb [6| 
CDF 0, 

D0 i. 



(4) 



In the absence of a world average, we use the most precise 
measurement from LHCb. For Atos the average of the 
CDF 3 and LHCb [13 measurements is 



Arris = ruH - rriL = (17.731 ± 0.045) ps~ 



(5) 



One should naturally ask if there are any constraints 
on the mixing parameters, beyond the obvious one: that 
the mass and width eigenvalues of the heavy and light 
mass eigenstates, tuh^l and Th.Li must be positive. (We 
use the notation customary in B physics, but the results 
apply equally for K'^ and D'^ mixing as well.) The time 
evolution of the flavor eigenstates is 
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(6) 



where M and F are 2x2 Hermitian matrices, and CPT 
invariance implies Mn = M22 and Fn = F22. The phys- 
ical states are the eigenvectors of the Hamiltonian, 

\BH,L)=p\B")Tq\B'), (7) 



where we chose |pp -I- \q\'^ — 1. CP violation in mixing 
occurs if the mass and CP eigenstates do not coincide. 



{Bh\Bl) 



b|2 - |g|2 _ 1 _ |5/p|2 



bP + kP i + k/p|' 
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The solution for the mixing parameters satisfies 
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(8) 



(9) 



p2 2A/i2 - iFi2 ' 
and from this and Eq. (j8]) it follows that (see, e.g., \Xm,) 
|2Afi2| |Fi2 



6 < min 



|ri2 



|2Mi 



(10) 



The measurable CP asymmetry in semileptonic (or any 
"flavor-specific" ) decay can be expressed as 



A 



SL 



1 - \q/p\' 
1 + \q/p\' 



25 



6^ 



Im (F12/M12) 
-|Fi2|V(4|Mi2 



(11) 
Thus, in the smaU S limit, Asl = 2(5-1- 0{S^). 

In the IF12/M12I ^ 1 limit, which applies model inde- 
pendently for the i?° ^ systems. 



Am 
AF 



2|Mi2|[l + 0(|Fi2/Afi2n], 

2 IF12I cos[arg(-Fi2/Mi2)] [l + 0(|Fi2/Mi2n] , 

Im (Fi2/Afi2) [1 -I- 0(|Fi2/A/i2n] . (12) 

In this limit, Eq. (|9]) implies that q/p is a pure phase 
to a good approximation, determined by Afi2, which has 
good sensitivity to NP. However, if IF12/M12I = 0(1), 
relevant for K'^ and D'^ mesons, then q/p depends on 
both F12 and Afi2 and the sensitivity to NP in Afi2 (and 
in argA/12) is diluted T2]. In that case Eq. ([T^ does not 
hold, but ^sL = 2(5 is a good approximation even in the 
D^ system, where the current bound is \S\ < 0.2 [J| . 

An additional constraint, the unitarity bound |13lll4|. 
stems from the time-evolution of the normalization of any 
linear combination of |i3°) and \B'~') being determined 
entirely by the F matrix. As discussed below, this con- 
strains the eigenvalues of F to be positive definite inde- 
pendent of the physical eigenvalues, or equivalently 



S'< 



ThT 



Hl^ L 



(mH-mi)2 + (Fff+Fi)V4 1 + x^ 



(13) 



Here we define, using T = {Th + Tl)/2, 



f 



r, 



y 



H 



2r 



(14) 



where x is positive by definition, wliile y £ (—1, +!)■ 

One may ask if other constraints exist purely from 
consistency considerations. As we show in a sepa- 
rate paper [15|, no additional limit on S exists, either 
physical or mathematical, without some knowledge of 
the Hamiltonian. In p articular, the invalidity of the 
bounds claimed in [la, ll3| can be made apparent by 
introducing new short-distance physics, which reduces 
I cos[arg(ri2/Mi2)]| while it leaves \Mi2\ and ri2 un- 
changed. 

However, this does not preclude the presence of new 
relations from appropriate inclusions of theoretical pre- 
dictions from models of the underlying interactions. In 
this paper, we derive a relation between the mixing pa- 
rameters of a meson system and only the magnitude |ri2 1 , 
which is in tension with the D0 measurement in Eq. ([T} . 
It has been known that the data in Eqs. ([T|) - ([3]) is 
not only in tension with the SM, but — assuming the 
SM calculation of ri2 — also with all models in which 
NP enters only through M12 [la Si- This is because 
Eqs. ([To]) and ([T2|) imply S < |ri2|/Am. Our result goes 
beyond this, because it makes optimal use of data on AF 
without theoretical assumptions, and indicates a larger 
tension independent of the nature of typical new physics. 



II. UNITARITY WITH THEORY INPUT 



As mentioned above, Eq. ([131 was first derived in 
Refs. [IJ, [ij]. We show that a stronger bound on S can 
be obtained using additional input from theoretical cal- 
culations. An analogy to the derivation of Ref. [Tj] will 
be particularly useful in deriving our results. 

We define the complex quantities 

a, = y^2^,{f,\H\B) , a, = y^2^, {f,\n\B) , (15) 

with Pi denoting the phase space density for final state 
fi. If we treat a^ and a^ as vectors in a complex iV- 
dimcnsional vector space, then taking the standard inner 
product on complex vector spaces, and using the optical 
theorem Il3ll. amounts to the relations 



a,- Qi 



Til 



To 



ri2 



(16) 



where CPT fixes Fn = F22 = F. Applying the Cauchv- 
Schwarz inequality to the vectors a^ and a^ implies [l3J 

|ri2|<Fn. (17) 

This is equivalent to the statement that the eigenvalues of 
the F matrix must be positive (in addition to Th,l > 0). 
To see that this is also equivalent to the unitarity 
bound of Eq. (IT3[) . we use Eq. ([7|) to define new vectors 
aHi and a^i analogously, such that 

^i = TT ^"-Hi + o-Li) , at = -— (ttLi - am) ■ (18) 

2p 2g 



For these newly defined vectors we can derive, in a similar 
manner as for Eq. ([r5[). the relations 



'-Hi 



aHi = Fh , 



''Li 



aLi — Fi 



a^i aLi 



-i{mH — TTiL + iV) 6 . 



(19) 



Substituting Eq. ([THj) into Eq. ^, using Eq. ([Hj) and 
the \q/p\'^ = (1 — (5)7(1 + 5) identity, we obtain Fn and 
F12 in terms of Mh.l, ^h,l and d. The unitarity bound 
in Eq. (J13l) then arises from using the new expressions 
for Fii and F12 in Eq. (IT71) . 

The preceding derivation made no assumptions on the 
actual values of the matrix elements appearing in the 
problem. For the kaon system, for which this bound was 
originally derived, this was a necessity due to the dom- 
inance of long-distance physics in the result. For B^s 
mesons, the large mass scale mi, ^ Aqcd allows Fn and 
F12 to be calculated in an operator product expansion, 
and at leading order IF12/F11I = ©[(AQCD/^b)^ (ISt^^)]; 
where IGtt^ occurs due to a one-loop difference between 
the two calculations. (In the Bd system there is an ad- 
ditional CKM suppression.) Thus it makes sense to con- 
sider some theory input, and we define 



2/12 = IF12I/F. 



(20) 



Using this relation between these matrix elements and 
proceeding with the same steps as above, we obtain, in- 
stead of an inequality as in Eq. (|17p , 



S'^ 



yh 



yh 



|Fi2p-(AF)V4 
|Fi2|2 + (Am)2 



(21) 



This equation follows from the solution of the eigenvalue 
problem, and was previously derived in Ref. 11] with 
the resulting bound on S noted. ^ (It also follows from 
Eqs. (9) and (12) in [19].) For fixed x and y, S^ is mono- 
tonic in 2/12, so an upper bound on yi2 gives an upper 
bound on \S\. For j/12 < 1 the usual unitarity bound in 
Eq. ([T3[) is recovered. 

Equation ([^T|) can also be obtained from a scaling ar- 
gument: As S only depends on mixing parameters, it is 
independent of the value of F. One can then scale F by 
j/12, which cannot affect S but changes x — > x/yi2 and 
y -^ y/yi2- Eq. ([2T|) follows then from this argument 
and Eq. ([T3|). The derivations above make the physical 
origin of this relation clear. Even if CPT is violated, the 
scaling argument, and therefore Eq. ([2Tt holds, although 
Eq. ([IHI) is modified. This applies for \S\'^, as CPT vio- 
lation allows 6 to be complex. 

Even if a precise calculation of F12 is not possible or one 
assigns a very conservative uncertainty to it, an upper 
bound on yi2 implies an upper bound on \S\, which is 
stronger than that in Eq. ([T3|). For small values of j/12, 
as in the Bd system, this bound can be much stronger. 



^ We were unaware of this in vl of this paper, and we thank Luis 
Lavoura and Joao Silva for bringing this to our attention. 
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FIG. 1: Upper bounds on ^|l as a function of AF^, setting 
AFd = 0. The darker [lighter] shaded region is allowed using 
the 1(7 [2cr] range of the theory calculation of IFj^j"!- The 
pair of horizontal [vertical] lines show the la range of the 
measured J^slI from D0 [AFs from LHCb]. The other curves 
are described in the text. 



III. COMPARING DATA AND THEORY 

We can compare the absolute value of Ag^ measured 
by D0, with the result implied by the relation above. At 
present Eq. (|2ip only provides an upper bound on |(5|, as 
the uncertainties of ri2 and AF allow the numerators for 
both Bd and Bg to vanish. Denoting this upper bound 
by S'^l^ and using the weight factors from Eq. ([2]), 

|4lI < (1.188±0.044)<,,+ (0.812±0.044)<5i;_. (22) 



Since Eq. (^J) only bounds |(5|, the bound on Ag^ is 
not sensitive to possible cancellations between Ag^ and 
v4gL (cf., the opposite signs of Ag^ in the SM, although 
I^slI ^ I^slD- ^^ Ato^.s are precisely known, we plot 
the bound as a function of AF^g. If LHCb measures 



A; 



SL 



Ai^m 



then the above bound with modified co- 



efficients apply for that measurement. 

In Fig. [1] we set AF^ = 0, which gives the most con- 
servative bound. The darker shaded region shows the 
upper bound on A\^^ using the \a ranges for |r]^2^| in 
the SM [U, 2|rf2| = (0.087 ± 0.021) ps-i and 2|Ff2| = 
(2.74 ± 0.51) X 10-3 ps-i. The dashed [dotted] curve 
shows the impact of using the 2cr region for Ff2 [rf2], 
and the lighter shaded region includes both 2cr regions. 
The vertical boundaries of the shaded regions arise be- 
cause jAFsl > 2|Ff2| is unphysical. A tension between 
the ^gL measurement and the bound is visible, indepen- 
dent of the discrepancy between the D0 result and the 
global fit to the latest available experimental data ^2] . 

We derived not an absolute bound in the fashion of 
the unitarity bound but a relation between calculable 
and measured quantities. It is thus worth clarifying the 
relationship of our result to the stated 3.9 cr disagreement 
of ^II with the SM reported in [l[. 

The SM prediction of ^sl uses the calculation of F12, 



and IF12I also enters our bound; thus, the discrepan- 
cies are correlated. Although the calculation of IF12I 
and Im(Fi2) both rely on the same operator product ex- 
pansion and perturbation theory, the existence of large 
cancellations in Im(Fi2) may lead one to think that the 
uncertainties could be larger in its SM calculation than 
what is tractable in the behavior of its next-to-leading 
order calculation ^23, 24]. The sensitivity of F12 to new 
physics is generally weaker than that of M12 (see [2^ [20] 
for other options). Thus, it is interesting to determine 5 
from Eq. (PT|) . besides its direct calculation. 

Of course, independent of the D0 measurement of ^II, 
we can also compare the bound implied by our relation 
to the individual best bounds on the semileptonic asym- 
metries in Eq. ^. To this end, in Fig. [2] we plot Ag^ 



vs. AFrf (and similarly for Bs) allowed by Eq. (|21l) and 
the \a and 2cr ranges of the SM calculation of ]Fi2l [2 If . 
Here, there have been no discrepancies claimed between 
the theory predictions and measurements, but our rela- 
tion allows us to place a bound tighter than the current 
experimental constraints which is more robust than the 
purely theoretical SM calculation as outlined above. 

Using AFs from LHCb in Eq. (JH), and neglecting AF^, 
we find the 2cr level bounds, 



]4lI<7.4x10- 



1A^lI< 4.2x10" 



(23) 



While this bound on ^g^ may seem to disagree with 
Fig. [21 note that in the plot the uncertainties of Ff 2 and 
AFs are not combined. Propagating the uncertainties, 
IF22P — (Ars)^/4 is negative at the la level, an unphys- 
ical result, hence the 2cr bounds in Eq. (P5)) . This bound 
on AgL is better than the current best bound in Eq. ([3]) 
by more than a factor of 3, while that for Ag^ is compa- 
rable. (However, in the case of Bd this is driven primarily 
by the uncertainty in the lifetime difference. If a non-zero 
value of AFrf were observed, a better bound could be de- 
rived.) It is worth emphasizing that this implication goes 
in both directions, given that an observation of ^g^ ^ 
may happen before that of AF^ ^ 0. Due to Eq. (|2ip . 
as soon as one of the two is measured to be nonzero, the 
other is constrained to be significantly smaller at worst 
and given a definite prediction at best. 



IV. CONCLUSIONS 

We provided a physical derivation of Eq. ((2T|) for neu- 
tral meson oscillation parameters, especially relevant for 
the i?° ^ systems, which allows incorporating theoreti- 
cal input on lri2l without any approximation, and with 
or without CPT conservation. This input is typically 
insensitive to the nature of NP and avoids the largest 
uncertainties of the direct theoretical calculation of CP 
violation in mixing. The application to the two neutral 
B systems, taking into account the recent LHCb mea- 
surement [0], leads to bounds on the semileptonic CP 
asymmetries of both system. These bounds are in ten- 
sion with the D0 measurement of ^|l, while providing 
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FIG. 2: Left plot: the region allowed by Eq. ((2T]) in the ^sl - AFd plane. The SM calculation of \V'\;^\ at Icr [2cr] gives the 
darker [lighter] shaded region. Right plot: same for j4gL — AFs; the straight lines show the Ict range of the LHCb result for AFs. 



a bound on the individual asymmetries at comparable or 
better levels than the current experimental bounds. Ad- 
ditionally, once an unambiguous determination of AgL 
or AF is made, we can use it to constrain the other ob- 
servable. Refinements of the A^^ measurements are an 
important part of the future B physics program [23, |2^ 
to search for new physics at both LHCb and the e+e" 
_B factories. Future bounds will in particular be helpful 
to constrain the individual measurements of AgL against 
the SM as well as consistency checks. 
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